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Abstract 

£> ■ We study the evolution of radiating and viscous fluid spheres as- 

. suming an additional homothetic symmetry on the spherically sim- 

\ metric space-time. We match a very simple solution to the symmetry 

OO ' equations with the exterior one (Vaidya). We then obtain a system of 

. two ordinary differential equations which rule the dynamics, and find 

\ a self-similar collapse which is shear-free and with a barotropic equa- 

*q | tion of state. Considering a huge set of initial self-similar dynamics 

O -1 states, we work out a model with an acceptable physical behavior. 

U 

1 Introduction 

X 

Often many authors assume spherical symmetry and perfect fluid approxi- 
mation to face the problem of self-gravitating and collapsing distributions of 
matter. Also, they use extensively progressive waves or similarity solutions 
(see [1, 2] and references therein). If the fluid is perfect the only equation 
of state compatible with self-similar fluids is the barotropic one [2]. The 
present paper concerns in part with the validness of the barotropic equation 
of state for a viscous and radiating fluid sphere. 
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In general, there are two types of self-similar space-times depending on 
whether they are invariant or not under scale transformations. Scale-free 
self-similar solutions are the similarity solutions of type one and the resulting 
space-time admits homothetic Killing vectors. Type two similarity solutions 
are not invariant under the simple scaling group [3]-[6]. The self-similar 
symmetry has been reported to characterize these two types of self-similar 
space-times [7]. 

Spherically symmetry and homothetic space-times show naked singular- 
ities. Assumption of similarity rather than spherical symmetry is crucial 
in determining the nature of the singularity in any gravitationally collaps- 
ing configuration [8, 9]. So far, self-similar space-times have been studied 
mainly in cosmological contexts [10]-[15]. 

Considering that the perfect fluid approximation is likely to fail, at least 
in some stages of stellar collapse, in this paper we study radiating and viscous 
fluid spheres. Specifically, we have been concerned with the radiative shear 
viscosity and its effect on the gravitational collapse [16]— [18]. We do not 
consider here the temperature profiles to determine which processes can take 
place during the collapse. For this purpose, transport equations have been 
proposed to avoid pathological behaviors (see for instance [19] and references 
therein). The motivation of this work was a recent study of radiating and 
dissipative spheres [20]. We assume an additional symmetry (homothetic 
motion) within the viscous fluid sphere without heat flow in the streaming 
out limit. 

The organization of this paper is the following. Section 2 shows the field 
equations, the junction conditions and the surface equations. In section 3 we 
write the homothetic motion equations in a convenient form. We propose a 
very simple solution in section 4 to work out some models. Finally, in section 
5, we draw conclusions. 

2 Dynamics and matching 

2.1 Field equations 

To write the Einstein field equations we use the line element in Schwarzs 
child-like coordinates 

ds 2 = e u dt 2 - e x dr 2 - r 2 (d6 2 + sin 2 9d0 2 ) . (1) 
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where v = v(t, r) and A = \(t, r), with (t, r, 8, <fi) = (0, 1, 2, 3). 

In order to get physical input we introduce the Minkowski coordinates 
(r,x,y,z) by [21] 



dr = e u / 2 dt, dx = e x ^ 2 dr, dy = rd6, dz = r sin ddcf), 



(2) 



In these expressions v and A are constants, because they have only local 
values. 

Next we assume that, for an observer moving relative to these coordinates 
with velocity uj in the radial (x) direction, the space contains 

• a viscous fluid of density p, pressure p, effective bulk pressure p^ and 
effective shear pressure p v , and 

• unpolarized radiation of energy density e. 

For this moving observer, the covariant energy tensor in Minkowski coor- 
dinates is thus 



/ p + e -e \ 

— e p + e — p^ — 2p v 

P~ Pc+Pv 

V P-Pt+P v / 



(3) 



Note that from (2) the velocity of matter in the Schwarzschild coordinates 



is 



dt 



(4) 



Now, by means of a Lorentz boost and defining p = p — p^, p r = p — 2p v , 
p t = p+Pr, and e = e(l + cj)/(l — uj) we write the field equations in relativistic 
units (G — c — 1) as follows: 



p + p r uj^ 
1-uj 2 



+ e = 



8nr 



1 

e 

r 



- - A, 



(5) 



Pr 



puj 



+ e 



87rr 



+ v. 



(6) 



3 



1 2 

6Zn r 

e- v [2A >tt + A >4 (A >4 -i/ >t )]} (7) 

ip + Pr) T ^ + e=-^-e~^ (8) 
1 — uj z Sirr 

where the comma (,) represents partial differentiation with respect to the 
indicated coordinate. Equations (5)-(8) are formally the same as for an 
anisotropic fluid in the streaming out approximation. 

At this point, for the sake of completeness, we write the effective viscous 
pressures in terms of the bulk viscosity (, the volume expansion 0, the shear 
viscosity rj and the scalar shear a [18] 

pc = ce (9) 



Pv = (10) 



where 
= 



(1-W 2 )V2 

and 



,,/nfXt ^00 + \ \lo( v r u r 2oj\ 



in) 



" = ^{z- — VT^s) < 12 > 

We have four field equations for six physical variables (p, p, e, u, ( and 
77) and two geometrical variables [y and A). Obviously, we require additional 
assumptions to handle the problem consistently. First, however, we discuss 
the matching with the exterior solution and the surface equations that govern 
the dynamics. 

2.2 Junction conditions 

We describe the exterior space-time by the Vaidya metric 

dg2 = / _ 2^M(X) \ ^ + 2dudR _ r2 ^ 2 + gin2 Qd ^ (i3) 
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where u is a time-like coordinate so that u = constant represents, asymp- 
totically, null cones open to the future and R is a null coordinate {qrr = 0). 
The relationship between the coordinates (t,r,#,0) and (u,R,9,(j)) is 



u 



*-r-2Mln(^-l), J R = r (14) 



The exterior and interior solutions are separated by the surface r = a(t). 
To match both regions on this surface we require the Darmois junction con- 
ditions. Thus, demanding the continuity of the first fundamental form, we 
obtain 

= 1 - ^ (15) 

and 

Va = -A a (16) 

From now on the subscript a indicates that the quantity is evaluated at 
the surface. Now, instead of writing the junction conditions as usual, we 
demand the continuity of the first fundamental form and the continuity of the 
independent components of the energy-momentum flow. This last condition 
guarantees absence of singular behaviors on the surface. It is easy to check 
that [18, 22] 

Pa = PC a + 2Pr la (17) 

which expresses the discontinuity of the radial pressure in presence of viscous 
processes. 

2.3 Surface equations 

To write the surface equations we introduce the mass function m by means 

of 

e -A(r,t) = i_2m( r ,i)/r (18) 
Substituting (18) into (5) and (8) we obtain, after some arrangements, 



dm , 
— — = -4irr z 
dt 



^fp r + e(l - - 2m/r) 1 /V /2 
dt 



(19) 



This equation shows the energetics across the moving boundary of the fluid 
sphere. Evaluating (19) at the surface and using the boundary condition (17) 
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(which is equivalent to p Ta = 0), the energy loss is given by 



m a = -47ia 2 e a (l-2m a /a)(l-u a ) (20) 

Hereafter overdot indicates d/dt. The evolution of the boundary is governed 
by equation (4) evaluated at the surface 

a = (1 - 2m a /a)uj a (21) 

Scaling the total mass m a , the radius a and the time-like coordinate by the 
initial mass m a {t = 0) = m a (0), 

A = a/m a (0), M = m a /m a (0), t/m a (0) -> t 

and defining 

n = u a (23) 
E = 47ra 2 e a (l-n) (24) 
the surface equations can be written as 

A = Ftt (25) 

F=j[(l-F)n + 2E] (26) 

Equations (25) and (26) are general within spherical symmetry. We need a 
third surface equation to specify the dynamics completely for any set of initial 
conditions and a given luminosity profile E(t). For this purpose we can use 
equation (7) or appeal to the conservation equation Tf" = evaluated at the 
surface. But we follow here another route, that is, we assume that the space- 
time admits a one-parameter group of homothetic motion generated by a 
homothetic Killing vector orthogonal to the four-velocity. These assumptions 
introduce some restrictions on the surface equations as is shown in the next 
section. 
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3 Homothetic motion 



We assume that the spherically symmetric space-time within the fluid 
admits a one-parameter group of homothetic motions. In general, a global 
vector field £ on the manifold is called homothetic if £^g = 2ng holds on a 
local chart, where n is a constant on the manifold, and £ denotes the Lie 
derivative operator. If n ^ 0, £ is called proper homothetic and it can always 
be scaled so to have n — 1; if n — the £ is a Killing vector on the manifold 
[23]-[25]. So, after a constant rescaling we write 

= 2g (27) 
where the vector field £ has the general form 

Z = A(r,t)dt + r(r,t)d r (28) 
After simple manipulations we obtain from (27) 

T = r (29) 

A, r = (30) 

Am tt + Tm ;r = m (31) 

Ai/ t + IV r + 2A = 2 (32) 

We further assume that the four-velocity is orthogonal to the orbit of the 
group 

u = - e ("-*>/ 2 (33) 
r 

Thus we obtain a connection between the time-like component of the homo- 
thetic Killing vector and the surface variables, 

Mt) = f (34) 

Now, expanding v near the surface, using (15), (16), (34), and evaluating at 
r = a the equations (5), (8), (31) and (32), after straightforward manipula- 
tions we find the surface equation 

Q= [ -^^(3F-1-2E) (35) 

From now on we disregard the bulk effective pressure to promote algebraic 
consistence. 
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4 Modeling 



In order to work out models we define the self-similar variables 

X = ™ (36) 
r 

and 

Y = -e u ' 2 (37) 
r 

Thus, equations (31) and (32) read 

AX tt + rX r = (38) 

and 

AY t + rY r = (39) 

In general these equations have solutions X = X(q) and Y = Y(q) , where q 
is 

^ = re -Jdt/A ( 40 ) 
We propose the specific solutions 

X = dq k (41) 

and 

Y = C 2 q l (42) 

where Ci, C2, k and / are constants. 

Solutions (41) and (42) are restricted by (15) and (16). Therefore the 
geometrical variables are 

( r \ fc+i 
a) (43) 

/ r \ 2(/+l) 

,-F(!) (44) 
In order to get the unique luminosity 

E=]-[F(k + 2l + 3)-(k + l)] (45) 
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we use equations (5), (6), (43) and (44) together with the boundary condi- 
tions (15), (16) and (17) to find 

n = 2( t fl)t/-l) (46 » 

where 

Z = [F 2 (5k 2 + Akl + lOJfe + 4/ 2 + 12/ + 9) 



-2F(k + 1) (5k + 21 + 3) + 5k 2 + 6k + 1] * (47) 

Note that "+" in the numerator of (46) represents the collapsing solution 
and "— " an expanding one. We consider here only Q + situations. 

Now, combining equations (35) and (46) we obtain an equation f(F, k, I) = 
0, which is too lengthy to present here, but which permits us to model dif- 
ferent situations. The first one is the shear-free and self-similar collapse for 
which k = I = 0, m/a ~ 0.3096 (m and a are linear with time) and p = 
at any space-time point. The second possibility appears upon solving for 
/ = l(F(t = 0),k) and includes the previous case. For k ^ we obtain 
shearing models but the homothetic symmetry is broken for t > 0. 

We work out a "tricky" third scenario by "forgetting" the origin of pa- 
rameter I, proposing that it depends on time in a very special way. If we 
imagine N initial self-similar states which represents the history of the col- 
lapsing surface, the symmetry equations (31) and (32) are satisfied at every 
point of the space-time without taking into account the variation with time 
of I. Therefore, we integrate numerically only equations (25) and (26), with 
(45), (46) and with I = l(t). Here we use standard Runge-Kutta (fourth 
order) methods and the initial conditions 

A(0) = 3.255; F(0) w 0.3856 

Once the boundary evolution and its energetics are determined, we use (43) 
[or (18)] and (44) to calculate the physical variables from the field equations. 
Figures (l)-(4) sketch the ratio p/p, dr/dt, e and rj, respectively, for k = 
(2 )10~ 3 . These self-similar spheres do not have a barotropic equation of state 
[figure (1)]. All shells evolve with decreasing collapsing velocities [figure 2)]. 
This behavior seems to be connected with the absorption of energy shown in 
figure (3) in the late stage. Shear viscosity increases initially with collapse 
but later decreases with time on any shell. 
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5 Conclusions 



We have assumed an additional symmetry to the space-time, homothetic 
motion, to generate non-static and simple solutions. These solutions were 
matched with the Vaidya one. We found that self-similar spheres with a 
barotropic equation of state (p = 0) are shear-free, this result is in complete 
accord with theoretical expectation [2], [26] [27]. Other self-similar scenarios 
are possible as well if we assume the evolution of the surface as a huge set 
of initial self-similar states. The shear viscosity profiles obtained in this 
work coincide qualitatively surprisingly well with others calculated in a more 
realistic framework [19]. 
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Time 



Figure 1: p/p as a function of time, for different values of r/a: 0.1 (uppermost 
curve), 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0 (lowermost curve). 
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Figure 2: dr/dt as a function of time, for different values oir/a: 0.1 (upper- 
most curve), 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0 (lowermost curve). 
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Figure 3: e as a function of time, for different values of r/a: 0.1 (initially up- 
permost curve), 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0 (initially lowermost 
curve). 
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4.0 




Time 



Figure 4: rj as a function of time, for different values of r/a: 0.1 (uppermost 
curve), 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0 (lowermost curve). 
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